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We present the first simultaneous radio, X-ray, ultraviolet, and optical spectroscopic observations of 
the M8.5 dwarf TVLM 513-46546, with a duration of 9 hours. These observations are part of a program 
to study the origin of magnetic activity in ultracool dwarfs, and its impact on chromospheric and coronal 
emission. Here we detect steady quiescent radio emission superposed with multiple short-duration, highly 
polarized flares; there is no evidence for periodic bursts previously reported for this object, indicating 
their transient nature. We also detect soft X-ray emission, with Lx/^boi ~ 10~ 49 , the faintest to date 
for any object later than M5, and a possible weak X-ray flare. TVLM 513-46546 continues the trend of 
severe violation of the radio/X-ray correlation in ultracool dwarfs, by nearly 4 orders of magnitude. From 
the optical spectroscopy we find that the Balmer line luminosity exceeds the X-ray luminosity by a factor 
of a few, suggesting that, unlike in early M dwarfs, chromospheric heating may not be due to coronal 
X-ray emission. More importantly, we detect sinusoidal Ha and H/3 equivalent width light curves with a 
period of 2 hr, matching the rotation period of TVLM 513-46546. This is the first known example of such 
Balmer line behavior, which points to a co-rotating chromospheric hot spot or an extended magnetic 
structure, with a covering fraction of about 50%. This feature may be transitory based on the apparent 
decline in light curve peak during the four observed maxima. From the radio data we infer a large scale 
and steady magnetic field of ~ 10 2 G, in good agreement with the value required for confinement of the 
X-ray emitting plasma. A large scale field is also required by the sinusoidal Balmer line emission. The 
radio flares, on the other hand, are produced in a component of the field with a strength of ~ 3 kG and 
a likely multi-polar configuration. The overall lack of correlation between the various activity indicators 
suggests that the short duration radio flares do not have a strong influence on the chromosphere and 
corona, and that the chromospheric emission may not be the result of coronal heating. 

Subject headings: radio continuum:stars — stars:activity — stars:low-mass, brown dwarfs — 
stars:magnetic fields 



1. INTRODUCTION 

In recent years it has become evident that low mass stars 
and brown dwarfs (spectral classes late-M and L) are capa- 
ble of producing unanticipated levels of magnetic activity, 
manifested primarily in their strong quiescent and flaring 
radio emission (Berger et al. 2001; Berger 2002; Berger 
et al. 2005; Burgasser & Putman 2005; Berger 2006; Os- 
ten et al. 2006b; Antonova et al. 2007; Phan-Bao et al. 
2007; Hallinan et al. 2007; Audard et al. 2007). Other tra- 
ditional activity indicators such as chromospheric Balmer 
line emission and coronal X-ray emission, however, appear 
to steeply decline from peak levels in early and mid M 



dwarfs (Pallavicini et al. 1981; Vilhu & Walter 1987; Gizis 
et al. 2000; West et al. 2004), leading to divergent trends of 
magnetic activity at the bottom of the main sequence. In 
both Ha and X-rays there is also a clear transition from 
persistent emission to a small number of flaring objects 
with duty cycles of a few percent (Reid et al. 1999; Gizis 
et al. 2000; Rutledge et al. 2000; Liebert et al. 2003; West 
et al. 2004), as well as a breakdown of the rotation-activity 
relation (Basri & Marcy 1995; Mohanty & Basri 2003) that 
is clearly seen in early M dwarfs (Rosner et al. 1985; Flem- 
ing et al. 1993; Mohanty et al. 2002; Pizzolato et al. 2003). 
The change and divergence in activity trends is most 
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clearly evident in the breakdown of the radio/X-ray cor- 
relation that holds for a large number of early-type stars 
and solar flares (Guedel & Benz 1993; Gudel et al. 1993a; 
Benz & Guedel 1994), and is attributed to flare heating 
of coronal plasma to X-ray temperatures (Neupert 1968; 
Guedel et al. 1996). While objects in the range M0-M6 
obey this correlation, several objects later than M7 ex- 
hibit radio emission that is several orders of magnitude 
brighter than expected (Berger et al. 2001; Berger 2002; 
Berger et al. 2005; Berger 2006). Similarly, in early M 
dwarfs there is an overall energy balance between X-ray 
and chromospheric emission, which has led to the idea of 
chromospheric heating by coronal X-rays (e.g., Cram 1982; 
Hawlcy et al. 1995). It is not known whether this mech- 
anism holds in ultracool dwarfs, primarily because of the 
decline in persistent activity. 

Theoretical work on magnetic dynamos in ultracool 
dwarfs also remains inconclusive. Studies of the a 2 dy- 
namo in fully convective stars suggest that a stratified 
and rotating turbulent medium can lead to the build-up of 
a non-axisymmetric and multi-polar field (e.g., Chabrier 
& Kiiker 2006; Dobler et al. 2006), but these models 
make several simplifications for computational purposes. 
It has also been argued that decreasing electrical conduc- 
tivity will impede the dissipation of any magnetic fields in 
the cool and increasingly neutral atmospheres of ultracool 
dwarfs (Mohanty et al. 2002). The existing radio detec- 
tions suggest that field dissipation may not be a problem, 
but the overall field configuration and the effect of neu- 
tral atmospheres in the presence of magnetic dissipation 
remain largely unexplored. 

As a result of the various conflicting trends and the 
transition from persistent to flaring emission, progress 
in our understanding of magnetic activity and dynamos 
in ultracool dwarfs requires simultaneous observations of 
the various activity bands. We have already under- 
taken such observations for the L3.5 brown dwarf 2MASS 
J00361617+1821104 in late 2002 and discovered periodic 
radio emission (P — 184 min), with no corresponding X- 
ray or Ha emission (Berger et al. 2005). These observa- 
tions indicated a magnetic field of <~ 200 G covering a sub- 
stantial fraction of the stellar surface, as well as the first 
direct confirmation that the radio/X-ray correlation is in- 
deed violated by orders of magnitude. Follow-up observa- 
tions showed that the field is stable on a ;> 3 yr timescale, 
much longer than the convective turnover time, pointing 
to a stable dynamo process. A recent simultaneous obser- 
vation of the L dwarf binary Kelu-1 resulted in an X-ray 
detection without corresponding radio emission (Audard 
et al. 2007), although the radio limits still allow for a vio- 
lation of the radio/X-ray correlation by up to <~ 2 x 10 3 . 

Here we exploit the powerful approach of simultaneous 
observations to investigate the magnetic activity in the 
M8.5 dwarf TVLM 513-46546, an object previously de- 
tected in the radio (Berger 2002, 2006; Osten et al. 2006b; 
Hallinan et al. 2006, 2007) and in Ha (Martin et al. 1994; 
Reid et al. 2002; Mohanty & Basri 2003). These observa- 
tions are the first in a scries that targets several objects in 
the sparsely-studied and critical spectral type range M7 to 

12 The VLA is operated by the National Radio Astronomy Obscrv 
crative agreement by Associated Universities, Inc. 



L3. In the case of TVLM 513-46546 we detect radio, X-ray, 
and Balmcr line emission, but no UV emission. The over- 
all behavior is complex and largely uncorrelated between 
the various emission bands. The long time and wavelength 
baselines provide unprecedented detail, including the first 
case to date of sinusoidal Balmcr line emission; the ob- 
served 2-hour period is in excellent agreement with the 
rotation of TVLM 513-46546. Using the various activity 
indicators we infer the properties of the magnetic field, 
corona, and chromosphere, and show that the underlying 
processes and field configuration likely differ from those in 
early M dwarfs. 

2. OBSERVATIONS 

We targeted the M8.5 dwarf TVLM 513-46546 due to its 
vicinity (d = 10.6 pc; Dahn et al. 2002) and known radio 
and Ha activity. The bolomctric luminosity of TVLM 513- 
46546 is Lboi ~ io -3 - 59 Lq, and its rotation velocity is 
vs'mi 60 km s^ 1 (Mohanty & Basri 2003). The non- 
detection of lithium, with a limit of 0.05 A, suggests that 
TVLM 513-46546 is most likely a very low mass star (Reid 
et al. 2002). Adaptive optics imaging of TVLM 513-46546 
revealed no companions with 5m <; 3 mag in the range 
0.1-15" (Close et al. 2003). 

TVLM 513-46546 was first detected in the radio during 
a 2 hr observation at 8.5 GHz in Sep. 2001, and exhib- 
ited both persistent {F v w 190 /iJy) and flaring emis- 
sion (Berger 2002), the latter with a peak brightness of 
1 mJy, a duration of 15 min, and circular polarization of 
r c w 66%. Subsequent observations from 1.4 to 8.5 GHz 
in Jan. 2004 revealed a similar level of persistent emission, 
F„(1.4) w 260 ,uJy, F v (4.9) w 280 ^Jy, and F„{8.5) w 230 
/iJy, with r c <; 15% (Osten et al. 2006b). Observations in 
Jan. 2005 revealed brighter emission, F„(4.9) w 405 /uJy, 
and F„(8.5) w 400 /iJy, as well as a claimed periodicity of 
about 2 hr (Hallinan et al. 2006). Finally, observations in 
May 2006 uncovered a series of flares with durations of a 
few minutes, ~ 100% circular polarization, and a period- 
icity of 1.96 hr (Hallinan et al. 2007). 

Previous detections of Ha emission reveal long term 
variability, with equivalent widths (EW) ranging from 1.7 
to 3.5 A, or L Ka /L hol w 10~ 5 (Martin et al. 1994; Reid 
et al. 2002; Mohanty & Basri 2003). 

Our simultaneous observations were obtained on 2007 
April 20 UT for a total of 8.8 hr in the radio (04:00-12:48 
UT), 8.9 hr in the X-rays (03:47-12:41 UT), and 7 hr in 
the optical (07:13-14:13 UT). Swift UV/optical telescope 
(UVOT) observations took place intermittently between 
03:58 UT and 12:10 UT, with a total on-source exposure 
time of 8036 s. 

2.1. Radio 

Very Large Array 12 observations were conducted at a 
frequency of 8.46 GHz in the standard continuum mode 
with 2 x 50 MHz contiguous bands. Scans of 295 s on 
source were interleaved with 50 s scans on the phase cali- 
brator J1513+236. The flux density scale was determined 
using the extragalactic source 3C48 (J0137+331). 

a facility of the National Science Foundation operated under coop- 
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The data were reduced and analyzed using the Astro- 
nomical Image Processing System (AIPS). The visibility 
data were inspected for quality, and noisy points were re- 
moved. To search for source variability, we constructed 
light curves using the following method. We removed 
all the bright field sources using the AIPS/IMAGR rou- 
tine to CLEAN the region around each source, and the 
AIPS/UVSUB routine to subtract the resulting source 
models from the visibility data. We then plotted the 
real part of the complex visibilities at the position 
of TVLM 513-46546 as a function of time using the 
AIPS/DFTPL routine. The subtraction of field sources is 
required since their sidelobes and the change in the shape 
of the synthesized beam during the observation result in 
flux variations over the map that may contaminate real 
variability or generate false variability. The resulting light 
curves are shown in Figures 1 and 2. 

2.2. X-Rays 

The observations were made with the Chandra/ ACIS-S3 
(backside-illuminated chip), with TVLM 513-46546 offset 
from the on-axis focal point by 15". A total of 29.76 ks 
were obtained. Data were analyzed using CIAO version 
3.3, and counts were extracted in a 1" radius circle cen- 
tered on the source position. We find a total of 8 counts 
in the 0.2 — 2 keV range, and 2 additional counts with 
kT 10 keV. Background counts were extracted from 
annuli centered on the source position, excluding other 
point sources detected in the observation. We find that 
2 background counts are expected within the source ex- 
traction aperture, likely corresponding to the two photons 
with kT w 10 keV. 

The source counts exhibit a narrow energy range with 
(kT) = 930 ± 250 eV, corresponding to a typical plasma 
temperature of 1.1 x 10 7 K. Using this temperature with a 
Raymond-Smith plasma model we find an energy conver- 
sion factor of 1 count = 3.4 x 10~ 12 erg cm~ 2 s _1 (0.2 - 2 
keV). Thus, the observed count rate of 2.69 x 10~ 4 s _1 
translates to a flux of 9.3 x 10~ 16 erg cuT 2 s" 1 (0.2 - 2 
keV), or a flux density of 1.4 nJy at kT = 1 keV. At the 
distance of TVLM 513-46546 the corresponding luminosity 
is Lx ~ 1.2 x 10 25 erg s _1 , or a ratio of Lx/^boi ~ 10~ 4,9 . 
This detection is at the same level as the quiescent emis- 
sion from the M8 dwarf VB 10 (Fleming et al. 2003), the 
faintest X-ray emitting late-M dwarf to date. 

We next find that of the 8 detected photons 4 arrive as 
pairs with separations of 217 and 31 s (Figure 1). The 
chance probabilities of such short time separations in a 
29.76 ks observation are 1.7 x 10~ 3 and 3.4 x 10~ 5 , respec- 
tively. It is thus possible that the second pair constitutes a 
flare. If true, the flare luminosity is 3.1x 10 24 erg cm~ 2 s _1 , 
or Lx/L\:, \ k, 10~ 5 ' 5 , the lowest luminosity flare detected 
from any late-M dwarf to date. The quiescent component 
would be correspondingly lower, Lx/L^oi ~ 10~ 50 . As we 
show below, the putative X-ray flare may coincide with the 
peak of the broadest radio flare. 

2.3. Optical Spectroscopy 

13 Observations were obtained as part of program GN-2007A-Q-60. 



We used 13 the Gemini Multi-Object Spectrograph 
(CMOS; Hook et al. 2004) mounted on the Gemini-North 
8-m telescope with the B600 grating set at a central wave- 
length of 5250 A, and with a 1" slit. A series of eighty 
300-s exposures were obtained with a readout time of 18 
s providing 94% efficiency. The individual exposures were 
reduced using the gemini package in IRAF (for bias sub- 
traction and flat-fielding), and rectification and sky sub- 
traction were performed using the method and software 
described in Kelson (2003). Wavelength calibration was 
performed using CuAr arc lamps and air-to-vacuum cor- 
rections were applied. The spectrum covers 3840-6680 A 
at a resolution of about 5 A. 

To measure the equivalent widths of the Ha and H/3 
emission lines we use continuum regions centered on 6551 
and 6572 A, and on 4854 and 4870 A, respectively. Sam- 
ple spectra in the low and high Balmer emission state are 
shown in Figure 3. The Ha light curve exhibits a clear 
sinusoidal behavior (Figure 1). 

2.4. Ultraviolet 

The data were obtained with the Swift /WOT in the 
UVW1 filter (Acs w 2510 A), as a series of 6 images with 
exposure times ranging from 560 to 1630 s (Figure 1). No 
source is detected at the position of TVLM 513-46546 in 
any of the individual exposures, or in the combined image 
with a total exposure time of 8036 s. We performed pho- 
tometry on the combined exposure using a circular aper- 
ture matched to the PSF of the UVW1 filter (2.2"), and 
found a 3ct limit of F A (UVW1) < 2.4 x 10~ 18 erg cuT 2 s" 1 
A" 1 , or a Vega magnitude of m(UVWl) > 23.0 mag. This 
limit corresponds to a ratio of UV to bolometric luminosity 
of \L\/ Lboi < lO- 3 ' 2 . 

3. MULTI- WAVELENGTH EMISSION PROPERTIES 

We observed TVLM 513-46546 across a wide wavelength 
range that traces activity in various layers of the outer 
atmosphere. The radio emission traces particle accel- 
eration by magnetic processes, and corresponds to gy- 
rosynchrotron radiation or coherent radiation (electron cy- 
clotron maser or plasma emission). The Balmer emission 
lines are thought to be collisionally excited in the chro- 
mosphere, and the X-ray thermal emission arises in the 
corona. 

3.1. Quiescent Emission 

Using the observed X-ray flux and temperature in the 
context of bremsstrahlung radiation we can estimate the 
coronal gas density and pressure. The emissivity is given 
by: 

r) v w 7.6 x 10~ 38 \n(T/v) n\ T~ 1/2 ergs" 1 cuT 3 Hz" 1 , 

where for TVLM 513-46546, T = 1.1 x 10 7 and v = 
2.2 x 10 17 Hz. Assuming a uniform coronal structure ex- 
tending w (0.1 - l)i?* w (0.7 - 7) x 10 9 cm above the 
photosphere, we find F„ t x ~ (0-3 — 6) x 10~ 52 n 2 erg s _1 
cm~ 2 Hz -1 . From the observed flux density of 1.4 nJy 
wc thus infer an electron density, n e w (0.5 — 2) x 10 10 



4 



cm~ 3 . We note that a coronal structure with a low vol- 
ume filling factor, /, would have a correspondingly higher 
density, n e oc / _1 ' 2 . Given this overall uncertainty, the 
inferred density is similar to what has been found for the 
M4.5 flare star AD Leo (flares: n e ~ 10 10 - 10 11 cm" 3 , 
quiescent: n e < 3 x 10 10 cm -3 ; van den Besselaar et al. 
2003) and is somewhat lower than in the M3.5 flare star 
EV Lac (n e ~ 10 13 cm" 3 ; Osten et al. 2006a). 

From the temperature and inferred density we find a 
coronal gas pressure of P ~ 10 — 40 dyne cm -2 . This does 
not take into account turbulent pressure, which at coronal 
temperatures may in fact be dominant (e.g., Osten et al. 
2006a. If the corona is confined by magnetic fields, the re- 
quired average field strength is thus B > \/8ttP ;> 15 — 30 
G. 

An alternative estimate of the coronal physical param- 
eters can be provided in the context of the loop model 
calculated by Rosner et al. (1978). In this model the 
loop temperature, density, and pressure are related by 
T w 1.4 x 10 3 (Pi) 1 / 3 K, where it is assumed that the loop 
apex temperature is similar to the mean coronal temper- 
ature. The loop length, £, is determined by the observed 
X-ray luminosity, I w 2.2 x 10 14 /T 7 / 2 L x x cm, where / is 
the filling factor of the loops. The simple relation between 
T, P, and I holds as long as P ~ const along the loop. 
This requires the loop length to be less than the pressure 
scale height, H w 7.5 x 10 2 T cm (for R = 0.1 R and 
M — 0.08 M Q ). From the observed parameters we find 
H w 8 x 10 9 cm and i w 8 x 10 13 / cm. The condition 
I < H thus requires a very small filling factor, / <; 10~ 4 . 
The pressure inferred from using t <~ H is P w 60 dyne 
cm -2 , and the magnetic field required for confinement is 
thus B J> 40 G, in good agreement with the value inferred 
above by assuming a uniform coronal structure. 

We now turn to the radio emission. The bremsstrahlung 
contribution in the radio band is F v w 0.03 /iJy, orders 
of magnitude lower than even the baseline quiescent emis- 
sion, which from several flarc-frcc regions of the radio light 
curve is found to be F„(8.5) = 208 ± 18 ^Jy. The 3cr 
limit on the fraction of circular polarization of the quies- 
cent component is r c < 25%. Both the flux and degree 
of circular polarization are similar to those measured in 
previous observations (§2), indicating that the quiescent 
component is stable on a multi-year timescale. 

Based on the brightness of the radio emission compared 
to the predicted thermal emission, and its long term sta- 
bility, we conclude that it is most likely due to gyrosyn- 
chrotron radiation. We follow the typical assumption that 
the mildly relativistic electrons, which produce the radio 
emission, follow a power law distribution, N("f) oc 7~ p for 
7 > 7m, with p ~ 3 typical for M and L dwarfs (Gudel 
et al. 1993b; Berger et al. 2005; Osten et al. 2006b). The 
gyrosynchrotron emission spectrum is determined by the 
size of the emission region (7?), the density of radiating 
electrons (n e ), and the magnetic field strength (B) ac- 
cording to (Dulk & Marsh 1982): 

r c = Q^xl0 1 - 93cose - 1 - 16cos2e {3xl0 3 /B)-°- 21 ~ - 37sine , (2) 

u m = 1.8 x 10 4 (sin^) 5 (n e i?) 23 B 077 Hz, (3) 
F v<m = 2.5 x 10~ 41 B 248 R 3 n e (sine)- 1 - 52 M Jy, (4) 



where 9 is the angle between the magnetic field and the 
line of sight. 

From previous observations of TVLM 513-46546 it ap- 
pears that the peak of the quiescent emission spectrum is 
v m k> 5 GHz (Osten et al. 2006b). Using the limit r c < 
25%, we infer a magnetic field strength, B <; 10—740 G, for 
6 = 20-80°. With this range we find R ~ (0.8-8) x 10*0 
cm, and n e ~ 100 — 1.3 x 10 10 cm~ 3 ; the latter range is 
for 9 = 80° to 20°. A comparison to the coronal density 
estimated from X-rays indicates 9 ~ 20 — 30°, and hence 
R ~ few x 10 10 cm and B <; few x 10 2 G. The inferred 
field strength is consistent with the value required for con- 
finement of the X-ray emitting coronal plasma as derived 
above. 

Finally, we turn to the observed Ha emission. The light 
curve exhibits clear sinusoidal behavior with a range of 
equivalent widths of 1.5 — 5.5 A (Figure 1). The sinu- 
soidal behavior indicates that the line is likely modulated 
by the rotation of TVLM 513-46546 rather than flares, and 
the emission is thus persistent in origin, at least on the 
timescale of our observation. Walkowicz et al. (2004) de- 
termined a multiplicative factor, \ = /x656o//boi, to con- 
vert Ha EW to L Ha /L hoh For TVLM 513-46546, the ob- 
served I — K = 4.3 mag indicates logx w —5.3, which 
matches the average value for spectral type M8.5 (Walkow- 
icz et al. 2004). Thus, for the full range of EWs we find 
log(LH a /-£'boi) ~ —4.6 to —5.1. This covers the typical 
range of Ha emission observed from M8.5 dwarfs (West 
et al. 2004) , as well as the range of EWs from past obser- 
vations of this source (§2). We return to the implications 
of the sinusoidal variations in §3.3. 

Since the X-ray and Ha emission appear to be persis- 
tent, we can gain insight into the physical processes in 
the outer atmosphere by examining the energy scale in 
each band. In early dMe stars there is evidence that the 
quiescent chromosphere is heated by downward directed 
coronal X-ray emission (Cram 1982). Nearly half of the 
absorbed energy will be radiatively lost in the Balmcr 
lines, primarily Ha (Cram 1982). Here wc find that the 
ratio L-ft a /Lx ~ 2 is in rough agreement with this overall 
model, taking into account the overall uncertainties in the 
X-ray and Ha luminosities. However, as noted by Cram 
(1982), any significant asymmetry in the coronal structure 
(e.g., extended loops) will substantially reduce the X-ray 
flux directed at the chromosphere. Since we find that the 
Ha line by itself already accounts for the maximal coronal 
heating rate, this model of chromospheric heating requires 
a uniform coronal structure. 

Alternatively, the marginal excess chromospheric lumi- 
nosity may indicate that the structure of the outer atmo- 
sphere in this late-M dwarf is markedly different than in 
early dMe stars. Potentially, the Ha emission is dominated 
by a non-chromospheric component. As we discuss below, 
in the context of the periodic emission, it is possible that 
the Ha emission region is a more extended "bubble" that 
co-rotates with the star. 

3.2. Radio Flares 

In addition to the persistent emission in the radio band 
we detect several distinct short duration (St ~ 2— 15 min) 
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flares, and a single broad brightening with a duration of 
about 1 hr (Figures 1 and 2). These flares range in peak 
flux density from 2 to 5.5 mJy. Our observations allow the 
detection of flares to a 5a sensitivity of about 3.5 mJy in a 
single 5 s integration. In Figure 2 we provide a zoom-in on 
the most distinct flares, including both the total intensity 
light curve and the circular polarization light curve. About 
half of the flares exhibit left circular polarization, while the 
other half exhibit right circular polarization, with overall 
fractions of ~ 50 — 100%. There is no apparent regular- 
ity in the duration, brightness, arrival time, or circular 
polarization of the flares. 

The apparently random sense of circular polarization 
(right- vs. left-handed) and flare arrival times suggest that 
the flares are produced in unrelated regions. This is con- 
trary to the model proposed by Hallinan et al. (2007) to 
explain their detection of periodic flares, with P = 1.96 hr 
well-matched to the rotation of TVLM 513-46546. These 
authors point to emission from a distinct region that is sta- 
ble over at least several rotation periods. The absence of 
a similar behavior in our data suggests that this stability 
is limited to < 1 yr. 

The variation in the sense of circular polarization from 
one flare to the next is also different from the uniform 
sense of polarization observed in radio flares from UV Cet 
(M5.5) and YZ CMi (M4.5). Flares in these two stars show 
consistent right-handed and left-handed polarization, re- 
spectively, which has been interpreted as a signature par- 
ticle acceleration in a toroidal dipole field (Kellett et al. 
2002). The behavior of radio flares from TVLM 513-46546 
in the current observation thus points to a multi-polar field 
configuration. 

While the quiescent radio emission is due to gyrosyn- 
chrotron radiation, the short durations and large frac- 
tion of circular polarization of the flares point to emis- 
sion by coherent processes. The electron cyclotron maser 
(ECM) process leads to emission at the fundamental cy- 
clotron frequency, v c = 2.8 x 10 6 B Hz, while plasma radi- 
ation is dominated by the fundamental plasma frequency, 
Up = 9 x 10 3 nj 2 Hz. In the ECM case the magnetic field 
strength inferred from the detection of flares at 8.46 GHz 
is B ss 3 kG, while in the case of plasma radiation we infer 
an electron density of n e ~ 9 x 10 11 cm~ 3 . The latter 
should be used as an upper limit in the ECM model. 

In the standard picture of the electron cyclotron maser 
(e.g., Melrose & Dulk 1982), electrons with a large pitch 
angle are reflected at the legs of the magnetic loop, while 
those with smaller angles precipitate into and heat the 
chromosphere. This process repeats, with the rise time 
of the emitted flare presumably reflecting the particle ac- 
celeration timescale. The flare decay reflects the time to 
precipitate out of the magnetic trap, and is thus related to 
the loss rate, r ~ [(Qz,/47r)i;/.£] _1 , where f^L is the solid 
angle of the loss cones, v is the electron velocity, and I is 
the transit length (Melrose & Dulk 1982). From our ob- 
servations at 8.5 GHz we infer v ~ c, and hence t ~ 30 s 
leads to, I ~ 10 12 (Ql/^k) cm. The actual size of the loop 
region is likely to be much smaller since the accelerated 
electrons traverse the field lines multiple times before the 
maser process ceases. Thus, whereas in the case of the 



periodic outbursts a stable, dipolar field configuration was 
inferred (Hallinan et al. 2007), here it is likely that the 
emission is arising from a tangled, multi-polar field. 

Finally, the ECM processes is supposed to provide coro- 
nal and chromosphcric heating as electrons precipitate out 
of the trap, and thus lead to increased Ha and X-ray emis- 
sion. From Figure 1 it is clear that this is not the case for 
TVLM 513-46546. None of the flares radio are accompa- 
nied by an increase in either the Ha or X-ray emission. 
This may provide further support to the proposed tangled 
field configuration, since it will lead to chromosphcric and 
coronal effects on a small scale that may go undetected 
when averaged over the stellar disk. 

3.3. Sinusoidal Ha Emission 

As noted above the Hen and H/3 emission lines exhibit 
sinusoidal variations with a periodicity of about 2 hours 
(Figure 1). This indicates that the likely origin of the 
Balmcr line modulation is the rotation of TVLM 513- 
46546. Indeed with w 7 x 10 9 cm, the 2-hour period 
requires a rotation velocity of v w 60 km s _1 , in perfect 
agreement with the observed vsmi w 60 km s _1 . If the 
Balmer lines are modulated by rotation, this indicates an 
inclination of the rotation axis relative to the line of sight 
of i w 90°. 

The likely connection to stellar rotation points to a large 
scale (/ ~ 50%) hot spot located on one hemisphere of 
TVLM 513-46546, or an emission "bubble" that may ex- 
tend to several times the stellar radius, and is occulted by 
the star for about half a rotation period. As noted above, 
the latter scenario may be relevant since the observed Ha 
luminosity appears to exceed the X-ray luminosity by a 
factor of a few, raising the possibility of an Ha emission 
region distinct from the chromosphere. The fact that the 
Balmer lines do not completely disappear during the light 
curve minima suggests that in both scenarios the solid an- 
gle subtended by the emission region is larger than the 
stellar disk, or alternatively, that the minima represent 
the baseline chromoshperic emission, while the rotation of 
an extended bubble or a hot spot into our line of sight 
produces the maxima. 

It also appears from the Ha light curve that the peak 
equivalent width decreases with time (Figure 1). The first 
peak has EW w 5.5 A followed by about 4.5 A for the two 
subsequent peaks, and about 3.8 A for the fourth peak. 
While our observations start during the decline of a pre- 
ceding cycle, an extrapolation back to the time of the ex- 
pected peak indicates an equivalent width of at least 6 A. 
The equivalent widths at the minima, however, appear to 
be more stable, with perhaps a slight decrease from 2 to 1.5 
A during our observation. These trends suggest that the 
hot spot or extended bubble may be transitory, possibly 
excited by an energetic event prior to the beginning of our 
observations. Extrapolating the observed trend, we find 
that the peak Ha equivalent width will match the baseline 
level of 1.5 A about 7 hrs after the end of our observa- 
tion. Future observations of TVLM 513-46546 are crucial 
for assessing whether this is in fact a transient feature. 

We further find no clear correspondence between the 
peaks of Ha emission and the radio flares. In fact, the 
opposite may be the case. The three Ha minima that 
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have simultaneous radio coverage appear to roughly coin- 
cide with radio flares. If we extrapolate the Ha light curve 
back in time, we find that a previous minimum also ap- 
pears to coincide with a radio flare (at about 05:40 UT). 
However, other radio flares, including the brightest one 
detected, do not coincide with observed or extrapolated 
Ha minima, suggesting that this anti-correlation between 
Ha and radio emission may be a result of small number 
statistics. 

The detection of a 2 hour period in two different emis- 
sion bands and at two different times with no clear cor- 
respondence indicates that the radio flares are not suffi- 
ciently energetic, or impact a large enough scale to influ- 
ence the chromosphere. This behavior also suggests that 
the magnetic field configuration of TVLM 513-46546 may 
be composed of a large-scale dipolar field, as well as a 
more compact and tangled component, which give rise to 
periodic signals at different times and wavebands. 

4. LACK OF RADIO /X-RAY CORRELATION 

TVLM 513-46546 is only the third ultracool dwarf to be 
observed simultaneously in the radio and X-rays, allowing 
a continued investigation of the radio/X-ray correlation in 
these objects. Coronally active stars up to spectral type 
M7, including the Sun, exhibit a tight correlation between 
their radio and X-ray emission (Guedel & Benz 1993; Bcnz 
& Guedel 1994). The persistent emission follows a linear 
trend, Lr w 3 x 10~ 16 Lx Hz -1 , which extends over 6 
orders of magnitude in Lr (Guedel & Benz 1993), while 
for flares the relation is Lr w 5.4 x 10~ 27 L x 37 Hz -1 over 
8 orders of magnitude in Lr (Bcnz & Guedel 1994); see 
Figure 4. 

From several previous (non-simultaneous) observations 
it has become clear that the above correlations break 
down in at least some ultracool dwarfs (Berger et al. 2001; 
Berger 2002; Berger et al. 2005; Berger 2006), with a clear 
transition occurring at spectral type M7 (Berger 2006). 
This conclusion was supported by our simultaneous ob- 
servations of the L3.5 dwarf 2MASS J00361617+1821104 
(Berger et al. 2005). Here we find continued evidence for 
this trend. Using the persistent radio and X-ray lumi- 
nosities we find Lr/Lx ~ 2.3 x 10~ 12 Hz^" 1 , a factor of 
about 7800 times larger than expected. This is similar to 
the level of excess radio emission observed in the previous 
late-M and L dwarfs. If we take into account the flaring 
emission, with a typical peak luminosity of Lr w 4 x 10 14 
erg cm~ 2 s _1 Hz -1 , and assume that the single X-ray pho- 
ton pair with Lx ~ 3 x 10 24 erg cm~ 2 s _1 is indeed a flare 
(§2.2), we find that for flares L R /L X « 1.3 x lO" 10 Hz" 1 , 
a factor of about 10 7 larger than expected. 

The radio/X-ray correlation has been interpreted in the 
context of the Neupert effect (Neupert 1968). In this sce- 
nario, the radio emission is produced when coronal mag- 
netic loops reconnect and create a current sheet along 
which ambient electrons are accelerated. The accelerated 
electrons in turn drive an outflow of hot plasma into the 
corona as they interact with and evaporate the underlying 
chromospheric material. The interaction of the outflow- 
ing plasma with the electrons produces X-ray emission via 
the bremsstrahlung process. This mechanism points to a 
causal connection between particle acceleration, which is 



the source of radio emission, and plasma heating, which 
results in X-ray emission. Thus, the X-ray thermal energy 
should simply be related by a constant of proportional- 
ity to the integrated radio flux. The breakdown in the 
radio/X-ray correlation, and the clear lack of correspon- 
dence between the radio flares observed on TVLM 513- 
46546 and its X-ray emission suggest that the heating 
of coronal material is generally inefficient in objects later 
than M7. This may be due to the short lifetime and small 
size of the flare emitting regions (i.e., in the case that 
they arise from highly tangled and multi-polar fields), or 
to lower efficiency of coronal heating so that the bulk of 
the coronal emission is at temperatures lower than kT <~ 1 
keV. 

5. DISCUSSION AND CONCLUSIONS 

We presented simultaneous radio, X-ray, UV, and 
optical spectroscopic observations of the M8.5 dwarf 
TVLM 513-46546 that probe magnetic activity and its in- 
fluence on the stellar chromosphere and corona. These 
observations are the first in a series of several objects that 
span the sparsely-studied spectral type range M7-L3, over 
which the magnetic activity appears to exhibit a change 
in behavior compared to early spectral types. We find 
that TVLM 513-46546 exhibits a wide range of both quies- 
cent and flaring activity that includes radio emission from 
large and small scale regions, coronal soft X-ray emission, 
and sinusoidal and periodic Balmer line emission from a 
chromospheric hot spot or an extended structure with a 
covering fraction of about 50%. 

Quantitatively, we find that the quiescent radio emission 
is produced on the scale of the entire stellar disk, with a 
magnetic field of ~ 10 2 G. From portions of the light curve 
that are free of flares we see no evidence for variability of 
this persistent component. Moreover, the similarity in flux 
level to observations carried over the past 6 years indicates 
that the magnetic field is stable on timescales significantly 
longer than the convective turnover time, T conv ~ 10 2 d 
(Chabrier & Baraffe 2000). This is similar to the situation 
we observed in 2MASS J00361617+1821104 (L3.5), with a 
periodic radio signal that was stable over at least 3 years 
(Berger et al. 2005). 

We also find from a comparison of the radio and X-ray 
data that the magnetic field axis is likely highly inclined 
relative to the rotation axis of TVLM 513-46546, which is 
inferred from the 2-hour period of the Ha emission to be 
about 90° (§3.3). This is an interesting result in the con- 
text of magnetic dynamo models of fully convective stars. 
Chabrier & Kiiker (2006) and Dobler et al. (2006) found 
that the a 2 dynamo, which relies on a stratified and rotat- 
ing turbulent medium, leads to a non-axisymmctric field 
with an overall configuration that lies in the equatorial 
plane. This seems to be supported by our observations. 

The X-ray emission requires a corona with T w 10 7 
K and a density of n e ~ 10 10 cm~ 3 , similar to those of 
early M dwarfs. The inferred coronal gas pressure re- 
quires a magnetic field strength of at least ~ 10 — 40 G 
for confinement, in good agreement with the radio-derived 
field strength. The energy input from the X-ray emitting 
corona is similar to, or somewhat smaller than, the radia- 
tive losses in the Balmer emission lines, indicating that 
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the chromosphere is at least partly heated by the over- 
lying corona. It is possible, however, that the somewhat 
elevated chromospheric luminosity is the result of an en- 
ergy input process that took place before the start of our 
observations. This latter possibility is supported by the 
apparently decreasing level of peak Ha flux during our ob- 
servation, and the nearly constant baseline level traced by 
the light curve minima. Indeed, the Ha luminosity during 
the minima is about half of the X-ray luminosity. 

In addition to the persistent emission, we detect a large 
number of radio flares with a range of peak fluxes, du- 
rations, and degrees of circular polarization. The over- 
all short durations of the flares and large degree of cir- 
cular polarization are indicative of coherent emission. In 
the context of the electron cyclotron maser mechanism, 
the inferred magnetic field is about 3 kG, similar to fields 
on the most active early M dwarfs (Saar & Linsky 1985; 
Johns-Krull & Valenti 1996). Similar flares have been de- 
tected in previous observations of TVLM 513-46546, but 
with a 2 hour periodicity that is absent in our data. The 
2 hour period was attributed to compact polar regions in 
a dipolar field rotating in and out of our line of sight (Hal- 
linan et al. 2007). The durations of the flares detected 
here, and their random arrival times and sense of circular 
polarization, point instead to a tangled and multi-polar 
field. Thus, the conditions required for coronal coherent 
radio flares exist on long timescales, but the change in 
behavior may signal a shift in the field configuration on 
<J 1 yr timescales. The inferred multi-polar nature of the 
field is again in good agreement with models of the a 2 dy- 
namo, which suggest that the bulk of the energy is in the 
quadrupolar and higher order components. 

Unlike in the radio flares, we do find clear Ha periodic- 
ity (P w 2 hr), with a sinusoidal light curve that reveals 
the presence of a chromospheric hot spot, or an extended 
bubble, with a covering fraction of about 50%. It is un- 
clear whether this emission region is stable over timescales 
longer than about 1 day, but the decrease in peak flux be- 
tween subsequent rotations may point to a transient nature 
that may be similar to the one now inferred in the radio 
band. The observed 2-hour period is well matched to the 
measured rotation velocity of TVLM 513-46546, and indi- 
cates a rotation axis inclination of about 90 deg. The ex- 
istence of such an extended structure provides additional 
support for a large-scale field that dominates the quiescent 
radio and X-ray emission. 

The general wisdom in the study of magnetic activity 
and its impact on the outer atmosphere is that the in- 
put of magnetic energy results in a series of related events 
that heat the corona and chromosphere and result in cor- 
related X-ray, radio, and optical line emission. This idea 
is supported by observations of the Sun, as well as various 
samples of early M dwarfs. The observations presented 
here show no clear evidence for any correlation between 
the various activity bands. In particular, the quiescent 
radio emission is over-luminous by nearly 4 orders of mag- 
nitude compared to predictions from the radio/X-ray cor- 
relation. Similarly, the radio flares do not appear to corre- 
late with the Ha variability. Finally, it is possible that the 
X-ray flux incident on the chromosphere is not sufficient 



to produce the observed Ha luminosity (particularly if we 
include the contribution from higher order Balmcr lines). 

Taking these various observations and inferences into 
account we therefore conclude that the observations of 
TVLM 513-46546 indicate that: 

• Ultracool dwarfs exhibit clear evidence for intense 
magnetic activity that is not diminished compared 
to early M dwarfs. This activity is manifested most 
clearly in the centimeter radio band, but also in X- 
rays and Balmer line emission. 

• Both a chromosphere and a corona appear to exist, 
with luminosities relative to the bolometric luminos- 
ity that are in line with other objects of similar spec- 
tral type, and are significantly lower (by ~ 2 orders 
of magnitude) compared to early M dwarfs. 

• The dissipation of magnetic energy leads to intense 
radio emission, but it does not have a clear effect 
on chromospheric and coronal emission, both tem- 
porally and in terms of overall luminosity. The 
radio/X-ray correlation is violated by about 4 orders 
of magnitude. 

• The decoupling of magnetic dissipation (as evi- 
denced by radio emission) from chromospheric and 
coronal emission may be due to changes in the struc- 
ture and scale of the magnetic field or the outer at- 
mosphere. This is possibly supported by the higher 
chromospheric luminosity compared to the coronal 
luminosity. 

• The presence of a steady large-scale magnetic field, 
as well as a multi-polar component, support current 
models of a 2 dynamos in fully convective stars. 

Simultaneous multi-wavelength observations of several 
additional ultracool dwarfs arc in progress. We expect 
that with this larger sample, and with the longer time 
baselines of our observations compared to typical studies, 
we can begin to address in detail the range of quiescent 
and variable activity, and the absence or presence of the 
correlations that appear to exist in the early M dwarfs. Ul- 
timately, these observations should reveal the energy and 
size scale of the magnetic field, and thus provide a detailed 
view of the magnetic generation process in fully convective 
stars. 
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Fig. 1.— Radio, Ha, UV, and X-ray of TVLM 513-46546. The arrival times of the X-ray photons are shown with arrows. 
Empty arrows correspond to likely background events with kT J> 10 keV; two are expected from the background count 
rate. Times of UV coverage are marked by gray squares; no UV emission is detected. The Ha emission is clearly sinusoidal 
and periodic, with P w 2 hr matching the rotation period of TVLM 513-46546, with an implied sini w 1. There is no clear 
correspondence between the various emission bands, with the possible exception of an X-ray photon pair that coincides 
with the broadest radio flare (at about 07:40 UT). 
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Fig. 2. — Zoom-in on individual radio flares. Total intensity (black) and circularly polarized flux (gray) are shown. The 
flares exhibit diverse behavior in terms of duration, amplitude, and fraction and sense of circular polarization. 
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Fig. 3. — Sample optical spectra of TVLM 513-46546 in the high and low Balmer emission line states. The high state 
spectrum has been offset upward for clarity. 




Fig. 4. — Radio vs. X-ray luminosity for stars exhibiting coronal activity. Data for late-M and L dwarfs are from Rutlcdge 
et al. (2000), Bcrger et al. (2001), Berger (2002), Berger et al. (2005), Burgasser & Putman (2005), Bcrger (2006), and 
Audard ct al. (2007), while other data are taken from Giidcl (2002) and references therein. Data for the Sun include 
impulsive and gradual flares, as well as microflarcs. The strong correlation between Lr and Lx is evident, but it begins 
to break down around spectral type M7 (see inset). 



